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The super-massive objects in galactic nuclei are thought to be the Kerr black holes predicted by 
General Relativity, although a definite proof of their actual nature is still lacking. The most massive 
objects in AGN (M ~ 10 9 Mq) seem to have a high radiative efficiency (r) ~ 0.4) and a moderate 
mass accretion rate (Lboi/^Edd ~ 0.3). The high radiative efficiency could suggest they are very 
rapidly-rotating black holes. The moderate luminosity could indicate that their accretion disk is 
geometrically thin. If so, these objects could be excellent candidates to test the Kerr black hole 
hypothesis. An accurate measurement of the radiative efficiency of an individual AGN may probe 
the geometry of the space-time around the black hole candidate with a precision comparable to the 
one achievable with future space-based gravitational-wave detectors like LISA. A robust evidence 
of the existence of a black hole candidate with r\ > 0.32 and accreting from a thin disk may be 
interpreted as an indication of new physics. For the time being, there are several issues to address 
before using AGN to test the Kerr paradigm, but the approach seems to be promising and capable 
of providing interesting results before the advent of gravitational wave astronomy. 
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I. INTRODUCTION 

Gravity has been tested and verified for distances in 
the range ~1 mm to ~1 pc (mainly within its Newtonian 
limit) and for weak gravitational fields [TJ [2J. The re- 
search is now moving to check the validity of the theory 
at cosmological scales, sub-millimeter distances, and for 
strong gravitational fields. One of the most intriguing 
predictions of General Relativity (GR) is that the col- 
lapsing matter produces singularities in the space-time. 
According to the weak cosmic censorship conjecture, sin- 
gularities of gravitational collapse must be hidden within 
black holes (BHs) [3]. In 4-dimensional GR, uncharged 
BHs are described by the Kerr solution, which is com- 
pletely specified by two parameters, the mass, M, and the 
spin angular momentum, J [3] . The condition for the ex- 
istence of the event horizon is a* < 1, where a* = \J/M 2 \ 
is the spin parameter^ When a* > 1, there is no horizon 
and the central singularity is naked, violating the weak 
cosmic censorship conjecture. 

Astronomers have discovered at least two classes of 
BH candidates (for a review, see e.g. Ref. [5]): stellar- 
mass objects in X-ray binary systems (M ~ 5 — 20 M@) 
and super-massive objects in galactic nuclei (M ~ 10 5 — 
10 9 Mq). The estimates of the masses of these objects 
are robust, because determined via dynamical measure- 
ments and without any assumption about the geometry 
of the space-time. The key-point is that the stellar-mass 
objects in X-ray binary systems are too heavy to be neu- 
tron or quark stars for any reasonable matter equation of 
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state [5], while the super-massive objects at the centers 
of galaxies are too heavy, compact, and old to be clusters 
of non-luminous bodies, as the cluster lifetime would be 
shorter than the age of these systems [7] . All these ob- 
jects are therefore thought to be the BHs predicted by 
GR, as they cannot be explained otherwise without in- 
troducing new physics. There are also some observations 
interpreted as an indirect evidence for the existence of the 
event horizon [S] (but see [§])■ On the contrary, there is 
no indication that the geometry around these objects is 
described by the Kerr metric. 

Testing the Kerr BH hypothesis is thus the next step to 
progress in this research field and several authors have in- 
deed suggested possible ways to do it using present and 
future data (for a review, see e.g. Ref. [ID])- A very 
promising approach is the detection of extreme mass ra- 
tio inspirals (EMRIs, i.e. systems consisting of a stellar- 
mass compact object orbiting a super-massive BH can- 
didate) with future space-based gravitational-wave an- 
tennas. Missions like LISA will be able to follow the 
stellar-mass compact object for millions of orbits around 
the central super-massive BH candidate, and therefore 
deviations from the Kerr geometry will lead to a phase 
difference in the gravitational waveforms that grows with 
the number of observed cycles [IT] . However, these data 
will not be available shortly, as the first mission will be at 
best in the early 2020s. The nature of BH candidates can 
also be tested by extending the methods currently used to 
estimate the spin of these objects, such as X-ray contin- 
uum |12j and Ka-iron measurements |13j , observations of 
quasi-periodic oscillations [H] , and measurements of the 
cosmic X-ray background [THl HE] . These methods can 
in principle be applied even with present data, provided 
that the systematic errors are properly understood. Fu- 
ture observations of the shadow of nearby super-massive 
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BH candidates are another exciting possibility to test the 
Kerr BH paradigm [17] . 

Previous studies have clearly pointed out that 
"rapidly-rotating" objects are the best candidates to test 
the Kerr BH hypothesis: if the object rotates fast, even 
a small deviation from the Kerr background can cause 
significant differences in the properties of the electromag- 
netic radiation emitted by the gas of the accretion disk 
and peculiar features, otherwise absent in the Kerr geom- 
etry, may show up jT5J Qj|] . The aim of this paper is to 
investigate the potentialities of the most massive BH can- 
didates in AGN, through the measurement of their radia- 
tive efficiency 77, defined by L^oi = 77M, where Lboi is the 
bolometric luminosity of the source and M is the mass 
accretion rate of the BH candidate. The estimate of the 
mean radiative efficiency of AGN through the Soltan's 
argument [2D] already suggests the presence of rapidly- 
rotating BHs [22 H3- Recently, Davis and Laor have 
proposed a way to measure the radiative efficiency of in- 
dividual AGN [23 . The authors found that the most 
massive BH candidates (with a mass M ~ 10 9 M@) 
would have a high radiative efficiency up to 7/ ~ 0.4, 
and a moderate mass accretion rate, iboi/^Edd ~ 0.3, 
where LecM is the Eddington luminosity of the source. 
The standard accretion disk model in a Kerr background 
would predict a high value of the spin parameter a* for 
these objects, extremely close to 1. At the same time, the 
moderate luminosity (in Eddington units) may indicate 
a thin accretion disk and the applicability of the stan- 
dard accretion disk model. If these estimates and these 
considerations are correct, the most massive compact ob- 
jects in AGN would be excellent candidates to test the 
Kerr paradigm. The sole measurement of rj can poten- 
tially constrain either a* and a deviation from the Kerr 
geometry. 

The paper is organized as follows. In Section [IT] I re- 
view the standard accretion disk model, its assumptions 
and properties, its effects on the evolution of the spin 
parameter of the central object, and its applicability. In 
Section 
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I consider an accretion disk in the Kerr back- 
ground, summarizing well-known results that should be 
expected if the BH candidates are the BHs of GR. In 
Sections [IV] and [Vj I discuss accretion disks respectively 
in the Johannsen-Psaltis (JP) [24] and in the Manko- 
Novikov (MN) 25j space-times. These are two metrics 
that can be conveniently used to describe a background 
deviating from the Kerr geometry. The nature of the two 
metrics is definitively different: the JP metric describes 
non-Kerr BHs in a putative alternative theory of grav- 
ity, while the MN one is an exact solution of the Ein- 
stein's vacuum equations and can describe the exterior 
gravitational field of generic compact objects. In both 
cases, the body is characterized by a mass, a spin angu- 
lar momentum, and an infinite number of "deformation 
parameters", even if here, for the sake of simplicity, I 
will consider only a single deformation parameter at a 
time. I will show that the two metrics present common 
features. In particular, a high radiative efficiency neces- 



sarily requires a very rigid compact object, much stiffer 
than a self-gravitating fluid with "normal" equations of 
state. The confirmation of the existence of individual 
AGN with high radiative efficiency (77 > 0.3) can poten- 
tially either be used to put strong constraints on the Kerr 
nature of astrophysical BH candidates and to discover 
new physics, as in the framework of the standard accre- 
tion disk model rj cannot exceed 0.32. In Section |VI| I 
discuss the findings of this paper in relation with current 
estimates of the radiative efficiency of AGN. The con- 
clusions and the issues that need to be addressed before 
using the most massive objects in AGN to really test the 



Kerr BH hypothesis are reported in Section VII In Ap- 
pendices [A] and [Bj the reader can find the non- vanishing 
metric coefficients respectively of the JP and of the MN 
metric. 



II. ACCRETION DISKS 
A. Novikov-Thorne model 

The Novikov-Thorne (NT) model is the standard 
model for accretion disks [55]. It describes geometrically 
thin and optically thick disks and it is the relativistic 
generalization of the Shakura-Sunyaev (SS) model [27] . 
The disk is thin in the sense that the disk opening angle 
is h = H/r -C 1, where H is the thickness of the disk 
at the radius r. Magnetic fields are ignored. In the Kerr 
background, there are four parameters (BH mass M, BH 
spin parameter a*, mass accretion rate M, and viscos- 
ity parameter a), but the model can be easily extended 
to any (quasi-)stationary, axisymmetric, and asymptoti- 
cally flat space-time. Accretion is possible because vis- 
cous magnetic/turbulent stresses and radiation transport 
energy and angular momentum outwards. The model 
assumes that the disk is on the equatorial plane and 
that the disk's gas moves on nearly geodesic circular or- 
bits. For long-term accretions, the disk is adjusted on 
the equatorial plane as a result of the Bardeen-Petterson 
effect j28j . That remains true even in non-Kerr back- 
grounds [29] . The assumption of nearly geodesic circular 
orbits requires that the radial pressure is negligible com- 
pared to the gravitational force of the BH. Heat advection 
is ignored (it scales as ~ h 2 ) and energy is radiated from 
the disk surface. 

The key-ingredient of the NT model is that the inner 
edge of the disk is at the innermost stable circular orbit 
(ISCO), where viscous stresses are assumed to vanish. 
When the gas's particles reach the ISCO, they quickly 
plunge into the BH, without emitting additional radia- 
tion. At first approximation, the total efficiency of the 
accretion process is 



?7tot = 1 — Eisco 



(2.1) 



where -Eisco is the specific energy of the gas at the ISCO 
radius and depends uniquely on the background geome- 
try. In general, the total power of the accretion process is 
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converted into radiation and kinetic energy of jet /wind 
outflows, so we can write rjtot = V + f}k- ?7 is the ra- 
diative efficiency and can be inferred from the bolomet- 
ric luminosity Lboi if the mass accretion rate is known: 
^boi = i]M. In this paper, I will assume 77k = 0, i.e. 
no gravitational energy of the gas is converted to kinetic 
energy of outflows. This is a conservative assumption for 
what follows. 

As a consequence of the accretion process, the BH spin 
parameter evolves. Since the gas particles arriving at 
the ISCO plunge quickly onto the central object, with- 
out emission of additional radiation, the BH changes its 
mass by SAI — iJisco^w and its spin angular momentum 
by SJ — LiscoSm, where I<isco is the specific angular 
momentum of the gas at the ISCO, while 8m is the gas 
rest-mass. The evolution of the spin parameter of the 
BH turns out to be governed by the following equation: 



da* 



1 ^isco 



dlnM M E] 



— 2a* 



ISCO 



(2.2) 



If the right hand side of Eq. (2.2 1 is positive, the accretion 



process spins the BH up. If it is negative, the BH is spun 
down. The equilibrium spin parameter a* q is reached 
when the right hand side of Eq. (2.2) vanishes and its 
value depends on the geometry of the space-ti me. As 
discussed in [T^ and explained bri efly a lso in Sec. IV the 
value of a* q we can infer from Eq. (|2.2| assuming that all 



the accreting matter has its angular momentum in the 
same direction provides the maximum value for the spin 
parameter of the super-massive BH candidates in galactic 
nuclei. This fact will be used to get the constraint a* < 



eq 
/ 1 



Actually, not all the radiation emitted by the disk can 
escape to infinity. A part of it leaves the disk, but it is 
then captured by the BH. Including this effect, Eqs. (2.1 ) 
and (2.2) become [30] 



dd: 



T) = 1 — Eisco — Ce 
1 ^isco + Cl 



dlnM M^isco + Ce 



— 2a* 



(2.3) 
(2.4) 



where Ce and Cl take into account the radiation captured 
by the BH and their expression can be found in (30] . 
There is also a part of radiation that is emitted by the 
disk and returns to the disk as a consequence of light de- 
flection (returning radiation) , interacting with the disk's 
particles and changing the flux profile emitted by the disk 
(see e.g. Ref. [31]). 



B. Validity of the model 

If we relax the assumption of vanishing stresses at the 
inner edge of the disk, the radiative efficiency increases. 
The bolometric luminosity of the disk becomes 



Lbo\ = n M = 

= Sisco^isco 



where gisco (gisco > 0) and Slisco are respectively the 
torque and the angular velocity of the gas at the ISCO. 
Another crucial point of the NT model is the a viscosity 
parameter: magnetic fields, which are thought to drive 
turbulence in disks, may not behave like a local scalar 
viscosity. If magnetohydrodynamics torques are present, 
the radiative efficiency 77 can even exceed 1, as the flow 
taps the spin energy of the BH [32 . In this case, it is 
not easy to recover -Eisco from the estimate of 77 and 
even the equilibrium spin parameter a* q is not given by 
Eq. J2~2l, but it turns out to be lower [33]. GRMHD sim- 



ulations show also that the gas's particles may not follow 
the geodesies of the space-time inside the ISCO [32] and 
that a significant emission of radiation from the plunging 
region is possible [35] . 

So, how good is the NT model to describe the ac- 
cretion disk around astrophysical BH candidates? For 
non-magnetized and weakly-magnetized disks, there is a 
common consensus that the NT model describes correctly 
thin disks, when the viscosity parameter is small, 

a <^ 1 |36J. A common criterion to select sources with 
thin disks is that L ho \/L^ AA does not exceed 0.3 [37] , 
In the case of magnetized disks (in Kerr background, as 
there are no simulations for other space-times), the issue 
is still open and controversial. The GRMHD simulations 
in Ref. [35] (see also [3J5]) show that the stress at the in- 
ner edge of the disk scales as h and the authors conclude 
that the NT model with vanishing stress boundary con- 
dition is recovered for h — > 0. The GRMHD simulations 
in [35] show instead large stress at the inner edge of the 
disk even when h 0; according to these authors, the 
NT model cannot describe magnetized disks, even when 
the disk is very thin. It is not clear if the disagreement 
between the two groups can be attributed to different 
configurations of the magnetic fields, different resolution 
of the simulations, or something else. 



III. KERR BLACK HOLES 

A Kerr BH is completely specified by two parameters: 
the mass, M, and the spin parameter, a*. In the NT 
model with 77k = 0, the radiative efficiency of the accre- 
tion process 77 is uniquely determined by the BH spin 
parameter, as a* sets the ISCO radius. There is a one- 
to-one correspondence between a* and 77, and 77 increases 
with increasing a^ At least in principle, an estimate of 
77 could be used to infer the spin parameter of the BH. In 
practice, assuming the validity of the NT model, we can 
deduce a lower bound for a*, as in general 77k 7^ 0. 

The radiative efficiency in Kerr background for differ- 
ent values of the BH spin parameter is shown in Tab. [I] 
and in the left panel of Fig. [I] If we use Eq. (2.1), for 



(1-£isco-Ce)M, (2.5) 



2 That is true only for Kerr BHs, i.e. when a„ < 1. It is not true 
in a Kerr background with arbitrary value of the spin parame- 
ter go]. 
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TABLE I. Kerr space-time. -Eisco, iisco, maximum radiative efficiency in the NT model, angular frequency at the ISCO 
radius, and ISCO radius in Boyer-Lindquist coordinates for different value of the spin parameter. * From Ref. |30) . without 
returning radiation, t From Ref. [31], including returning radiation. See text for details. 
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FIG. 1. Kerr space-time. Left panel: maximum radiative efficiency in the NT model, angular frequency at the ISCO radius, 
and ISCO radius (divided by 10) in Boyer-Lindquist coordinates as a function of the spin parameter. Units M = 1. Right 
panel: Evolution of the spin parameter for an initially non-rotating BH of mass Mq. See text for details. 



corotating accretion disks one finds that rj is in the range 
~ 0.057 (Schwarzschild BH, a* = 0) to ~ 0.423 (extreme 
Kerr BH, a* = 1). However, r) increases slowly with 
a* for low values of the spin parameters and it increases 
much faster when a* approaches 1. For instance, r\ > 0.30 
requires a* > 0.996, r\ > 0.34 requires a* > 0.999, and 
r\ > 0.40 demands a* extremely close to 1 (see Tab. [I] and 
Fig. [I}. 

In the Kerr background, we can analytically integrate 
Eq. f2~2j) and we obtain gT] 



2 M 
3M 
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.Ml 
M 2 



(for M/M < y/6) . 



a, = l (for M/M > V6) , (3.1) 
for an initially non-rotating BH with mass Mq- The 



equilibrium spin parameter is thus <z* q = 1, and it is 
reached after the object increased its mass by a factor 
V6 f=a 2.4. The evolution of the spin parameter as a 
function of M/M is shown in Fig. [l] right panel. 

For an astrophysical BH, the situation is slightly dif- 
ferent, as there is no realistic mechanism to spin up the 
object to a* too close to 1. The accretion process from a 
thin disk is still a very efficient mechanism, but we have 
to include the effect of the radiation emitted by the disk 
and captured by the BH, Eqs. (2.3) and (2.4 1. In this 
case, we find slightly lower values for r\ and ai H . However, 
as r] increases very quickly when a* approaches 1, even a 
tiny different in a* q causes a non-negligible difference in 
the maximum value of rj. As the radiation with angular 
momentum anti-parallel to the BH spin has larger cap- 
ture cross section, from Eqs. (2.4) one finds a* q < 1. For 
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instance, Thome found a° q = 0.9978 when the emission 
of the disk is isotropic, and a* q = 0.9983 when the emis- 
sion is limb-darkened [30]. The corresponding efficien- 



cies [Eq^ ( |2^3j )] are respectively 0.302 and 0.308. From 
Eq. pO] ), we would have found 0.3f8 (for a* = 0.9978) 
and 0.326 (for a* = 0.9983). The effect of the return- 
ing radiation, which changes a little bit the emission of 
the disk, introduces an even smaller correction |31j . The 
equilibrium spin parameter is now a* q = 0.9983 in the 
case of isotropic emission and at q — 0.9986 in the case of 
limb-darkened emission. The radiative efficiencies turn 
out to be resp ectiv ely 0.309 and 0.315 [0.326 and 0.331 
if we used Eq. (2.1)]. In conclusion, if the super-massive 



BH candidates are Kerr BHs, a realistic upper bound for 
i] should be around 0.32, as these objects can unlikely be 
spun up to a* higher than roughly 0.998. 

Tab. |T] shows also the angular frequency at the ISCO, 
^iSCOj and the value of the ISCO radius in Boyer- 
Lindquist coordinates, 7'iscch for the same values of the 
spin parameter. Even for rhsco and nsco there is a one- 
to-one correspondence (for a given mass M) with the spin 
parameter of the BH. In the case of corotating disks, the 
angular frequency at the ISCO increases with a*, from 
-0.068 M (a, =0) to 0.5 M (a, = 1). The ISCO radius 
decreases with a*, from 6 M (a* = 0) to M (a* = 1). 
^isco and nsco as a function of the spin parameter are 
shown (for corotating disk) in the left panel of Fig. [I] 



IV. JOHANNSEN-PSALTIS SPACE-TIMES 

The JP space-times have been proposed in [23] explic- 
itly to be used to test the Kerr geometry around astro- 
physical BH candidates. They describe BHs (at least 
for a* lower than a critical value) and they are an ex- 
tension of the Kerr solution, in the sense that here the 
compact objects are specified by the mass, the spin an- 
gular momentum, and an infinite number of deforma- 
tion parameters (i = 3, 4, 5, etc.) measuring devia- 
tions from the Kerr geometry. When all the deformation 
parameters vanish, one recovers exactly the Kerr back- 
ground. The explicit expression of the metric is reported 
in Appendix [X] The JP metric is not a solution of any 
known gravity theory, but it seems to be a very conve- 
nient framework to test the Kerr BH hypothesis. 

The NT model for thin accretion disks can be easily 
extended to non-Kerr backgrounds. The first important 
difference is the determination of the ISCO radius. In the 
Kerr background, circular orbits on the equatorial plane 
are always vertically stable and the ISCO is set by the or- 
bital stability along the radial direction. If the compact 
object is more oblate than a Kerr BH, that remains true. 
On the contrary, if the compact object is more prolate 
(including the case in which the object is still oblate, but 
simply less oblate than a Kerr BH with the same spin) 
both kinds of instabilities are possible. When the ISCO 
radius is marginally stable along the radial direction, as 
in Kerr, it is at the minimum of the energy of equatorial 



circular orbits. When instead the ISCO is marginally 
stable along the vertical direction, such a minimum does 
not exist in general. In the latter case, the curve of the 
energy of equatorial circular orbits around the ISCO is 
clearly steeper, causing a lower value of -Eisco with re- 
spect to the case of radially unstable ISCO with the same 
radius. The total efficiency r\ is therefore higher. 



As Eqs. (2.2 1 and (2.4) depend on the geometry of the 



space-time, the equilibrium spin parameter in these met- 
rics is determined by the deformation parameters ej. We 
can notice that in the case of the super-massive BH can- 
didates, a* q should be the maximum value for the spin 
parameter of these objects, independently of their actual 
nature [TBI [29] 142] . Indeed, for the super- massive BH 
candidates the initial value of a*, i.e. the one at the time 
of their birth, is thought to be completely irrelevant, as 
the mass of these objects has increased by several or- 
ders of magnitude from the original one and the spin has 
evolved accordingly. Long term accretion from a thin 
disk can efficiently spin the BH up to a* q (but it spins 
the BH down if a* > a* q ), while other processes (chaotic 
accretions, minor and major mergers) more likely spin 
the object down to a* — 0. So, the value of a* q we can 



infer from Eqs. (2.2) and (2.4) can be used to exclude the 



region a* > a* in the diagrams spin parameter vs de- 
formation parameter. This point will be crucial for what 
follows. 

Let us now start considering the properties of the JP 
space-time with deformation parameter £3 and £i = for 

i ^ 3. Fig. [2] shows the efficiency 77, as deduced from 

Eq. (2.1). In particular, we are interested in the most 



efficient systems and I report the curves on the de- 
plane with r) = 0.30, 0.32, 0.34, 0.36, 0.38, 0.40. The 
black solid curve marks the equilibrium spin parameter 
a* q . For £3 > 0, the compact object is more prolate than 
Kerr and a* q < 1. For £3 < 0, the object is more oblate 
and ai q > 1. It is thus clear that observations can po- 
tentially provide very strong constraints for £3 < 0, and 
much weaker bounds for £3 > 0. In particular, for £3 > 0, 
the maximum value of a* may be even significantly lower 
than 1, but still r\ can be very high. Interestingly, for 
£3 7^ 0, r] does not increases abruptly as a* approaches 
a* q and therefore the effect of the radiation emitted by 
the disk and captured by the BH may be neglected. In- 
deed, we can assume that the radiation captured by the 
BH decreases a* q by — 0.002, as in Kerr (actually the 
correction to a° q should be smaller than 0.002, as the 
ISCO radius is larger, see the right panel in Fig. k5). This 
is equivalent to replace the curve of a* q in Fig. 



the curve 



I with 

0.002. At this point, we have to decrease 
ij by — 0.01. Except for the Kerr space-time and for the 
space-times with small deviations from the Kerr geome- 
try, the two effects should not cause significant changes. 
In particular, we still find a large region with objects with 
77 > 0.32, which would be impossible for a Kerr BH in a 
realistic astrophysical context. 

Fig. [3] shows the contour plot of the angular fre- 
quency at the ISCO, ^isco, (left panel) and the one of 
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FIG. 2. JP space-time with deformation parameter £3 and ej = for i 7^ 3. Contour plots of the radiative efficiency 
7] = 1 — -Eisco: J? = 0.30 (red solid curve), 0.32 (blue dotted curve), 0.34 (red solid curve), 0.36 (blue dotted curve), 0.38 (red 
solid curve), 0.40 (blue dotted curve). The black solid curve is the equilibrium spin parameter a* q as inferred from Eq. (2.2 1 
and represents the maximum value for the spin parameter of the super-massive BH candidates. This means that the region 
with a* > aj q is not allowed. The right panel is simply the enlargement of the area inside the orange box in the left panel. 



0.4 




0.4 


0.3 




0.3 


0.2 

CQ 

CO 










0.2 

CO 

to 


0.1 








0.42 


0.1 







0.32\ 









-0.1 










-0.1 




0.92 0.94 


0.96 0.98 

a. 


1 


1.02 1.04 






0.92 0.94 0.96 0.98 1 1.02 1.04 

a. 



FIG. 3. JP space-time with deformation parameter £3 and ej = for i 7^ 3. Left panel: contour plot of the angular frequency at 
the ISCO radius: ^lisco/M = 0.32 (red solid curve), 0.34 (blue dotted curve), 0.36 (red solid curve), 0.38 (blue dotted curve), 
0.40 (red solid curve), 0.42 (blue dotted curve). Right panel: contour plot of the ISCO radius in Boyer-Lindquist coordinates: 
f^isco/M = 2.00 (red solid curve), 1.75 (blue dotted curve), 1.50 (red solid curve), 1.25 (blue dotted curve). The black solid 
curve is the equilibrium spin parameter u, q as inferred from Eq. (2.2 1. 



the ISCO radius in Boyer-Lindquist coordinates, nscOi 
(right panel). While accretion in JP space-time may have 
high efficiency, thus mimicking a Kerr BH with a* ex- 
tremely close to 1 (or even exceeding the efficiency of a 
Kerr BH with a* = 1), the same object would have f^isco 
and nsco more similar to the one expected for Kerr BHs 
with lower spin parameter, fiisco can be potentially in- 
ferred in several ways from observations (e.g. from the 
variability of the source) and that may be useful to break 
the degeneracy between a* and £3 in 77. 

Lastly, we can consider JP backgrounds with a dif- 
ferent deformation parameter. Interestingly, one finds 
essentially the same picture. In Fig. |4j I report the con- 
tour plots of r) for the JP space-time with deformation 



parameter e 4 and e, = for i ^ 4 (left panel) and for the 
one with deformation parameter £5 and = for i 7^ 5 
(right panel). For these two cases, even the contour plots 
of ^isco an d ?"isco ( n °t shown here) are fairly similar to 
the ones in Fig. [3] 



V. MANKO-NOVIKOV SPACE-TIMES 

The MN space-times are stationary, axisymmetric, and 
asymptotically flat exact solutions of the Einstein's vac- 
uum equations [25] , These metrics can describe the ex- 
terior gravitational field of a compact object with mass 
M, spin parameter a*, and arbitrary mass-multipole mo- 
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FIG. 4. As in Fig. [2] for the JP space-time with deformation parameter £4 and Cj = for i 7^ 4 (left panel) and for the one 
with deformation parameter £5 and £i = for i 7^ 5 (right panel). 




FIG. 5. As in Fig.gjfor the MN space-time with deformation parameter q. The compact object is oblate (prolate) when q > — 1 
(q < — 1), but it is more oblate (prolate) than a Kerr BH with the same value of the spin parameter for q > (q < 0). 




FIG. 6. As in Fig. [3] for the MN space-time with deformation parameter q. 
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ments, while the current-multipole moments are fixed by 
the former. The simplest non-Kerr object has three free 
parameters: mass M, spin parameter a*, and anomalous 
quadrupole moment q (see Appendix FBI) . q is defined by 



Q = (l + q)Q Kc 



(5.1) 



where Q is the mass-quadrupole moment of the compact 
object, while QKerr = —o^M 3 is the one of a Kerr BH 
with the same spin parameter a*. The case q = cor- 
responds to the Kerr metric, while for q > (q < 0) 
the object is more oblate (prolate) than a Kerr BH. The 
properties of this subclass of MN space-times with defor- 
mation parameter q were studied in [43| and |18j . When 
q < 0, for a quite restricted set of values of a* and q, the 
plunging region does not connect the IS CO to the com- 
pact object [18 . That occurs when inside the ISCO there 
is another region with stable circular orbits on the equa- 
torial plane and energy lower than Eisco ■ In these space- 
times, the gas's particles plunging from the ISCO get 
trapped between the ISCO and the compact object. As 
the gas need to radiate additional energy and angular mo- 
mentum in order to fall to the compact object, the value 
of the radiative efficiency turns out to be higher than 
the one inferred from Eq. (2.1) or Eq. (2.3). The actual 



value of r\ depends on the astrophysical processes. For 
instance, Ref. [IB] discussed a simple model in which the 
gas forms a thick disk inside the ISCO and r\ = 1 — E- m , 
where E ln is the specific energy at the inner edge of the 
thick disk, is up to a few percent higher than the one 
predicted by Eq. (2.1). 

Fig. [5] shows the radiative efficiency r\ as computed 
from Eq. ( |2.1[ ). The MN metric is written in prolate 
spheroidal coordinates, which require a* < 1; here the 
region a* > 1 is covered by the Malko-Mielke-Sanabria 
Gomez solution [33] which can be parametrized by the 
same three parameters of our subclass of the MN metric 
(i.e. M, a*, and q) and was shown to be very similar to 
the MN space-time in the common region of validity [32] • 
It is clear that there are strong analogies with the JP 
space-times. For objects more oblate than Kerr, a* q > 1 
and a high radiative efficiency is possible only for small 
deviations from the Kerr geometry. For objects more 
prolate than Kerr, a° q < 1 and when the ISCO radius 
is set by the orbital stability along the vertical direction 
7] can be high even for objects rotating slower. Fig. [6] 
shows the angular frequency at the ISCO, fiisco; an d 
the ISCO radius in Boyer-Lindquist coordinates, risco 
As in the JP background, the same object with high ra- 
diative efficiency has lower Qisco and larger nsco- The 
case of the MN space-time with q = and an anomalous 
mass-hexadecapole moment h shows similar properties in 
V, ^isco, and r Isco . 



VI. DISCUSSION 

The radiative efficiency of a source can be deduced 
from the formula rj — L\ m i/M. In general, however, 



it is not easy to get an estimate of 77, because a mea- 
surement of the mass accretion rate can be problematic. 
The Soltan's argument provides an elegant way to de- 
termine the mean radiative efficiency of AGN, fj, from 
the mean BH mass density in the contemporary Universe 
and the AGN luminosity per unit volume integrated over 
time |20j . There are several sources of uncertainty in 
the final result, but a conservative bound seems to be 
fj > 0.15 [21] . If the compact objects in AGN are Kerr 
BHs, this constraint requires a* > 0.89. In Ref. [22], the 
authors used a revised version of the Soltan's argument, 
in which some important assumptions are not necessary. 
They found a mean radiative efficiency fj sa 0.30 — 0.35 
in the redshift interval 0.4 < z < 2.1. 

Recently, it has been proposed a way to estimate the 
radiative efficiency of individual AGN [33] . The mass ac- 
cretion rate M can indeed be determined from the low 
frequency region of the thermal spectrum of the accre- 
tion disk of these objects, if the mass M is known. The 
standard accretion disk model cannot reproduce the ob- 
served spectrum of AGN. That is due to reprocessing, 
Comptonization in a corona, radiative transfer effects in 
the inner accretion disk, an so on. Here the key-point 
is that this method relies on the validity of the simple 
thin-disk model at relatively large radii, where these ef- 
fects are thought to be irrelevant. One can then estimate 
the total luminosity of the source, irrespective of its ex- 
act production mechanism, and infer the radiative effi- 
ciency. So, the redistribution of the radiation emitted in 
the accretion process does not affect the measurement of 
radiative efficiency; for more details, see Sections 2.3 and 
4.3 of Ref. [23] ■ The authors found a strong correlation 
of 77 with M, raising from r\ ~ 0.03 when M ~ 10 7 M 
and -Lboi/^Edd ~ 1 to 77 ~ 0.4 when M ~ 10 9 M G and 
-^bol/^Edd ~ 0.3. 

The most massive objects in AGN with M ~ 10 9 Mq 
may be excellent candidates to test the Kerr BH hypoth- 
esis. The standard criterion to consider thin disks is in- 
deed that iboi/^Edd does not exceed 0.3, corresponding 
to an opening angle h not larger than 0.05. In this case, 
the estimate of rj may provide a measurement of the spe- 
cific energy at the ISCO (or at least a lower bound, as 
we are neglecting r/^ and Ce)- The estimate of r\ can then 
be used to constrain possible deviations from the Kerr 
geometry and may have some advantages with respect to 
the most popular techniques to measure the BH spin pa- 
rameter, i.e. the continuum fitting method and the Ka 
iron line analysis. The continuum fitting method can be 
used only for stellar-mass BH candidates, for which there 
is no argument to constrain a* : the value of the spin pa- 
rameter of the BH candidates in X-ray binary systems 
should reflect the one at the time of the formation of the 
object and it is definitively impossible to predict its value 
relaxing the assumption that GR is the correct theory of 
gravity. The method of the Ka iron line can be used for 
the super-massive BH candidates in galactic nuclei (so 
we can impose a* < a* q ), but it necessarily relies on an 
ad hoc intrinsic surface brightness profile and there is no 
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way to determine it a priori. Although simulations have 
the potential to compute the surface brightness profile, 
the standard required to test the Kerr background is a 
much higher level of reliability than is likely to be met at 
any time in the foreseeable future. 

Let us now assume that future GRMHD simulations 
confirm the validity of the NT model for thin accretion 
disks and that future observations provide robust evi- 
dence that the most massive objects in AGN have a very 
high radiative efficiency and a moderate mass accretion 
rate. If we exclude the possibility of objects more pro- 
late than Kerr, which may be difficult to explain theoreti- 
calljj^J from the right panel of Fig. [5] we can see that q can 
be constrained at the level of ~ 1CP 3 . That roughly cor- 
responds to test the mass-quadrupole moment of these 



objects with a precision ~ 10 , see Eq. (5.1). For a 
self-gravitating fluid with reasonable equations of state, 
q is much larger, i.e. the compact object becomes sig- 
nificantly more oblate than a Kerr BH when the spin 
parameter increases. For instance, a neutron star should 
have q > 1 [45]. A similar bound may be obtained by 
observing the gravitational waves emitted by an EMRI 
with future space-based gravitational-wave detectors (see 
the third reference in [11]). even if a direct comparison 
is not possible, as here we are probing the geometry of 
the space-time very close to the compact objects, where 
higher orders multipole moments are not really negligible, 
while gravitational wave detectors will study the space- 
time at larger distances, and therefore they will be sen- 
sitive to the mass-quadrupole moment only. If we do not 
exclude a priori the existence of compact objects more 
prolate than a Kerr BH, the constraint is much weaker: 
for 77 > 0.30, deviations from the Kerr mass-quadrupole 
moments can still be up to 20%. Let us now consider a 
more exotic situation: the NT model works and observa- 
tions find unambiguously r\ > 0.32. In absence of torque 
at the inner edge of the disk, it would be difficult to 
explain this outcome in the Kerr background, and new 
physics may be invoked. Let us notice, however, that 
a non-Kerr background can explain high 77, say up to 
77 ~ 0.4, but higher values may not be natural, as we still 
have a* < a* q . 



3 The true problem is not the existence of compact objects less 
oblate than a Kerr BH with the same spin parameter, but that 
the objects we are talking about have the ISCO marginally sta- 
ble along the vertical direction. In the case of the JP metric, 
these objects are BHs with two disconnected event horizons, one 
above and one below the equatorial plane | 19 |. These objects 
more prolate than Kerr and high radiative efficiency have there- 
fore "two centers of attraction", one above and one below the 
equatorial plane. It is not clear at all if such a system can exist 
and be stable. 



VII. CONCLUSIONS 

There is some evidence that the most massive BH can- 
didates in AGN have a high radiative efficiency and a 
thin accretion disk. In this case, they could be excellent 
candidates to test GR in the strong field regime and, in 
particular, the Kerr BH paradigm. A robust observation 
of a high radiative efficiency together with the confirma- 
tion of the validity of the standard accretion disk model 
for the accretion process onto these objects would con- 
strain possible deviations from the Kerr geometry of the 
space-time around astrophysical BH candidates. If wc 
restrict our attention only to objects more oblate than a 
Kerr BH with the same spin parameter, a radiative effi- 
ciency 77 > 0.30 requires that the quadrupole moment of 
the compact object deviates not more than ~ 10~ 3 with 
respect to the one of a Kerr BH. If we allow for the ex- 
istence of compact objects more prolate than a Kerr BH 
(whose existence, however, may be questionable, see the 
footnote 3), the constraints is much weaker, at the level 
of 20%. On the other hand, the observation of an accret- 
ing BH candidate with 77 > 0.32 might indicate that the 
super-massive BH candidates in galactic nuclei are not 
the Kerr BH of GR, as there is no realistic astrophysical 
mechanism capable to spin up them to a* > 0.998. 

For the time being, there are at least two main issues 
to address before using the measurement of the radiative 
efficiency of individual AGN to test GR: 

1. We do not really know if the NT model can de- 
scribe thin accretion disks around astrophysical BH 
candidates. If it can or it is easy to estimate the 
necessary corrections, the measurement of 77 could 
provide an interesting way to probe the geometry of 
the space-time around BH candidates. If it cannot, 
it seems to be impossible to test GR (with this or 
other approaches) from the properties of the radi- 
ation emitted by the gas in the accretion disk and 
we should wait for the advent of gravitational wave 
astronomy. 

2. We need robust and more precise measurements of 
77. At present, there are several sources with 7/ > 
0.32 at high masses, but we cannot really exclude 
they actually have 77 < 0.32 [23] • There are some 
sources of uncertainty in the final estimate of the 
radiative efficiency and the most important one is 
the determination of the mass of the BH candidate. 
Moreover, the results presented in |23| are based on 
a preliminary study and it is necessary to further 
investigate and verify the validity of the method 
proposed by these authors. 
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Appendix A: JP metric 

The JP metric is not a solution in any known grav- 
ity theory. It is a simple parametrization to describe 
the space-time around non-Kerr BHs and was specifically 
proposed to test the Kerr BH hypothesis [24] . The metric 
was obtained by starting from a deformed Schwarzschild 
solution and then by applying a Newman-Janis trans- 
formation. The non-zero metric coefficients in Boyer- 
Lindquist coordinates are: 



9tt 



2Mr\ , 

i-^H (l + A) 



2aMr sin 2 



= sm 



2 i 2 

r + a 



(1 + h) , 

2a 2 Mr sin 2 



a 2 (p 2 + 2Mr) sin 4 



where 



9rr 
969 



P 

A 



P 2 (l + h) 



2 a ' 



A + a 2 h sin 2 9 

2 



r 2 + a 2 cos 2 I 



(Al) 



E 

k=0 



2Mr + a 2 , 
Mr 

£2k + 



"£2fc+l 



M 2 

„2 



(A2) 



The metric has an infinite number of free parameters 
and the Kerr solution is recovered when all these param- 
eters are set to zero. However, in order to recover the 



correct Newtonian limit we have to impose eo = e i = 0, 
while 62 is constrained at the level of 10~ 4 from current 
tests in the Solar System [2"4"] . 



Appendix B: MN metric 

The MN metric is a stationary axisymmetric, and 
asymptotically flat exact solution of the Einstein's 
vacuum equations with arbitrary mass-multipole mo- 
ments :25J. It does not describe the space-time around 
a BH and it has naked singularities and closed time-like 
curves at small radii. These pathological features should 
be either inside some sort of exotic compact object, whose 
exterior gravitational field would be described by the MN 
metric, or GR should break down close to them. The non- 
zero metric coefficients in prolate spheroidal coordinates 
are: 



9tt — —J , 

9t<)> = , 

94><t> = j( x2 - x ) i 1 - y 2 ) - f u 



9 XX 

9yy 



f 



fc 2 e 2 7 x 2 _ y 2 



(Bl) 



where 



/ 
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e 2 i'A(x 2 ~lY l (l-a 2 Y 2 . 



(B2) 
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A= (x — 1)(1 + ab) 2 - (1 - y 2 ){b - a) 2 
B=[x + l + (x- l)ab] 2 + [(1 + y)a + (1 - y)b] 2 , 

C =(x 2 - l)(l + ab)[b-a-y(a + b)} + (1 - y 2 ){b - a)[l + ab + x{\ - ab)] 
"+oo / n / n N 



a = —a cxp 



k=0 



b = a exp 



. n=l \ fc=0 /. 



(B3) 
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Here R = \J x 2 + y 2 — 1 and P n are the Legendre poly- 
nomials with argument xy/R, 

p - p (^L\ 

n ~ n \ R J ' 

The standard Boyer-Lindquist coordinates (r, 9) are re- 
lated to the prolate spheroidal coordinates (x, y) by 

r = kx + M , 
cos9 = y. (B5) 

The MN solution has an infinite number of free param- 
eters: k, which regulates the mass of the space-time; a, 
which regulates the spin; and a n (n = 1, ...,+oo) which 
regulates the mass-multipole moments, starting from the 
dipole a\, to the quadrupole a-i, etc. For a ^ and 
a n = 0, the MN solution reduces to the Kerr metric. For 
a = a n = 0, it reduces to the Schwarzschild solution. 
For a = and a n ^ 0, one obtains the static Weyl met- 
ric. Without loss of generality, we can put a± = to 
bring the massive object to the origin of the coordinate 
system. The simplest extension of the Kerr solution is 



thus the subclass of MN space-times with a n — for 
n =/= 2. Here there are three free parameters (fc, a, and 
o^) related to the mass M, the dimensionless spin pa- 
rameter a* = J/M 2 , and the dimensionless anomalous 
quadrupole moment q, defined by Q = — (1 + q)a 2 M 3 , 
with Q the mass-quadrupole moment of the object, by 
the relations 

VT^of-i 

a = ^ , 

a* 



a 2 = qa^-rp- . (B6) 

Note that q measures the deviation from the quadrupole 
moment of a Kerr BH. In particular, since Qkcit = 
— a 2 M 3 , the solution is oblate (prolate) for q > — 1 
(q < —1), but it is more oblate (prolate) than the Kerr 
one for q > (q < 0). When q — 0, the solution reduces 
to the Kerr metric, but when q ^ also the higher-order 
mass-multipole moments have a different value than in 
Kerr. 
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